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Magnetic exchange driven proximity effect at a magnetic insulator / topological insulator (MI/TI) 
interface provides a rich playground for novel phenomena as well as a way to realize low energy dis¬ 
sipation quantum devices. Here we report a dramatic enhancement of proximity exchange coupling 
in the MI/magnetic-TI EuS / Sb 2 -ccV x Te 3 hybrid heterostructure, where V doping is used to drive 
the TI (Sb 2 Tes) magnetic. We observe an artificial antiferromagnetic-like structure near the MI/TI 
interface, which may account for the enhanced proximity coupling. The interplay between the prox¬ 
imity effect and doping in a hybrid heterostructure provides insights into the engineering of magnetic 
ordering. 

PACS numbers: 61.05.fj, 75.25.-j, 75.30.Gw, 75.70.Cn. 


The time-reversal symmetry (TRS) breaking and sur¬ 
face bandgap opening of topological insulator (TI) is an 
essential step towards the observation of novel quantum 
phases and realization for Tl-based devices mm- In gen¬ 
eral, there are two approaches to break the TRS: transi¬ 
tional metal (TM) ion doping, where Cr or V are doped 
into the entire TI BHH], and magnetic proximity effect, 
where a magnetic insulator (MI) layer in proximity to TI 
induces exchange coupling |3j [6HE] . Doping TM impuri¬ 
ties into TI will introduce a perpendicular ferromagnetic 
(FM) anisotropy and provide a straightforward means 
to open up the bandgap of the TPs surface state, with 
profound influence to its electronic structure SUMS. 
In particular, quantum anomalous Hall effect (QAHE), 
where quantum Hall plateau and dissipationless chiral 
edge channels emerge at zero external magnetic field, 
has recently been realized in Cr-doped and V-doped TIs 
in nni nsHso]. Ideally, compared to the doping method, 
the proximity effect has a number of advantages, includ¬ 
ing spatially uniform magnetization, better controllabil¬ 
ity of surface state, and freedom from dopant-induced 
scattering that degrades TI properties, as well as pre¬ 
serving TI intrinsic crystalline structure, etc. eu na. 
However, due to the in-plane anisotropy and low Curie 
temperature of the common MI, such proximity effects 
are usually too weak to induce strong proximity mag¬ 
netism in TI. In fact, compared to magnetically doped TI 
which can induce as large as a 50meV surface bandgap 
[U, the EuS/TI system has only 7 meV gap opening due 


to the strongly localized Eu-f orbitals [23]. Therefore, the 
enhancement of proximity magnetism is highly desirable 
to make it a valuable approach as doping hence takes full 
advantage. 

In this Letter , we report significant enhancement of the 
proximity effect in MI EuS / magnetic-TI Sb 2 - cc V cc Te 3 
hybrid heterostructure. Using polarized neutron reflec- 
tometry (PNR), we inferred an increase of proximity 
magnetization per unit cell (u.c.) in TI, from 1.2/i^/u.c 
to 2.7/i#/u.c at x = 0.1 doping level. High-resolution 
transmission electron microscopy (HRTEM) image iden¬ 
tifies the TI/EuS interfacial sharpness and excludes false 
positive magnetism signal from interdiffused Eu ions into 
TI. Furthermore, the proximity effect enhancement is ac¬ 
companied by a decrease in the interfacial magnetization 
of EuS, resulting in an exotic antiferromagnetic (AF) like 
structure. The existence of the “artificial” AF structure 
between FM EuS and the FM Sb 2 -^V cc Te 3 with differ¬ 
ent anisotropies is consistent with magnetometry mea¬ 
surements which shows exchange bias. Such artificial AF 
ordering in this FM / FM hybrid heterostructure may 
shed lights for designing novel magnetic phases and de¬ 
vices. 

High-quality MI 6nm EuS / 15 quintuple layer (QL) 
magnetic TI Sb 2 - cc V ;E Te 3 hybrid heterostructures were 
grown by molecular beam epitaxy (MBE) under a base 
vacuum ^ 5 x 10“ 10 Torr, where magnetic TI thin films 
Sb 2 - cc V cc Te 3 were grown on clean, heat-treated sapphire 
(0001) substrates with V-dopants coevaporated during 
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FIG. 1. (color), (a) MI EuS / V-doped TI Sb 2 Te 3 hybrid het- 
erostructure. The arrows denote the spin direction. The V- 
doped TI layer has a perpendicular magnetic anisotropy, while 
the magnetic EuS has in-plane magnetic anisotropy. Such het¬ 
erostructure may create an exotic magnetic environment near 
the interface, as illustrated in (b). For a given Eu ion (red- 
arrow), it interacts with neighborhood intra-plane Eu (or¬ 
ange arrow) through Heisenberg interaction, inter-plane Eu 
ions (green arrow) through super-exchange interaction, spin- 
polarized states at TI surface and localized moments in TI 
(blue arrow). 



FIG. 2. (color), (a) The configuration of specular PNR. ki, kf 
and Q denote the incident wavevector, reflected wave vector 
and the momentum transfer, respectively, (b) The spin + and 
spin — PNR data R + and R~ for the EuS/Sbi.gVo.iTe 3 and 
EuS/Sb 2 Te 3 samples, at low (5mT) and high (0.7 T) in-plane 
guide fields. The fitting results are represented by solid lines 
and the curves are shifted for clarity, (c) The spin asymmetry 
for the reflectivity in b. (d) The same spin asymmetry, but as¬ 
suming one control sample EuS/Sb 2 Te 3 have exactly the same 
thickness as EuS/Sbi.gVo.iTe 3 . In this way, the spin asym¬ 
metry difference is dominated by only the magnetic structure 
instead of crystalline structure. At ~ 0.4 nm -1 , the differ¬ 
ence comes from the effect of external magnetic field, while 
at ~ l.Onm -1 the difference mainly comes from V-dopants. 


TI growth. The EuS (111) layer was deposited in situ 
over the TI film using electron beam source. To under¬ 
stand the interplay between proximity effect and TM 
doping, 6nm EuS/ 15 QL Sb 2 Te 3 , 15 QL Sb 2 - :E V ;E Te 3 
and 15 QL Sb 2 - cc V cc Te 3 samples were fabricated. 

The atomic configuration of the MI / magnetic TI het¬ 
erostructure is shown in Fig. [l] (a). The upper EuS 
has in-plane anisotropy [241127 : within xz-plane, while the 
lower TM doped TI has easy axis out-of-plane [lQl fl4l [20] 
along y-axis. The different anisotropy directions in the 
two layers, corroborated by a strong interfacial spin-orbit 
coupling, create a complex magnetic environment for the 
EuS near the interface (Fig. |T] b). The Heisenberg inter¬ 
action, superexchange interaction ^3 HI] , d — f coupling 
[29] and coupling with the TI’s spin texture may finally 
contribute to an overall augmentation of the proximity 
effect. 

The PNR experiments were carried out using PBR 
beamline at the NIST Center for Neutron Research, from 
which the in-plane magnetization is extracted. The ex¬ 
perimental setup is shown in Fig. |2ja), where the inci¬ 
dent spin-polarized neutron beam is reflected by the het¬ 
erostructure sample, while the spin non-flip reflectivity 

signals from both spin components (++ and-) were 

collected under external guide magnetic field. The PNR 


refinement is based on a multilayered TI / proximity /in¬ 
terfacial EuS / main EuS model [22] . To maximize the 
PNR information extraction, we did not compare the x 2 
with and without proximity effect due to limited sensitiv¬ 
ity. Instead, we presume the existence of the proximity 
coupling layer and optimize its magnitude. 

The spin non-flip reflectivity curves for the MI/ mag¬ 
netic TI sample EuS/Sbi.gVo.iTe 3 and control sample 
MI/TI EuS/Sb 2 Te 3 , at low (5mT) and high (0.7T) 
fields, are shown in Fig. [2] (b). The fitting and refinement 
of PNR is performed using the GenX program [30]. To 
directly infer the possible contribution of V-dopants, the 
corresponding spin asymmetries SA = are plot¬ 

ted (Fig. E i) for the raw data and thickness-adjusted 
data (Fig. |2i). In this way, the different features of the 
SA in Fig. 2] (d) are solely coming from the magnetic 
structure since the crystalline structure is taken to be 
identical. We see that at Q ~ 0.4 nm -1 , /iqH = 5mT 
SA for both samples with and without V-dopants over¬ 
lap each other, but distinct with the fi^H = 700 mT SA 
curves, indicating an effect from guide field within this 
Q range; while at Q ~ 1.0 nm -1 , a splitting of the SA 
curves for both samples at same guide field (eg. blue 
and red curves) is observed. This indicates the influence 
of the V-dopants to magnetic structure at high Q range 
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(spatially localized) even without performing fitting. 

Results of fitting the PNR data are shown in Fig. 3. 
The substrate lies in the region below Onm. Nuclear- 
SLD (NSLD, red curves) identifies the chemical compo¬ 
sitional contrast, where the NSLDs for each compound 
layer are correctly reproduced from PNR fitting (sapphire 
substrate 5.5 x IO -5 nm “ 2 , Sb 2 Te 3 1.8 x IO -4 nm “ 2 , 
EuS 1.5 x 10 -4 nm -2 and amorphous AI 2 O 3 capping layer 

4 x 10 - 4 nm -2 ). This further validates the fitting qual¬ 
ity. In Fig. |5Ja), without the EuS proximity layer, the V- 
dopants in the Sb 1 . 9 V 0 . 1 Te 3 sample contribute to no more 
than 0.2/is/u.c in-plane magnetization at /io H = 0.7 T, 
indicating a very strong perpendicular FM anisotropy in 
V-doped TI. This is consistent with the result in the in¬ 
set of Fig. [4] (a), and facilitates us in obtaining reliable 
PNR refinement by fixing the magnetization of the mag¬ 
netic TI layer. The magnetic-SLD (MSLD) (blue curves) 
at /io H = 5mT and 700mT in-plane guide fields at T = 

5 K are also plotted. In Fig. [ 3 ] (b), for the EuS/pure 
TI sample, we see a penetration of magnetization into 
TI. Unlike the EuS region where strong absorption LSD 
(ASLD) is always accompanied due to the Eu ions’ large 
neutron absorption, the magnetization into TI does not 
show any absorption (^ 14 — 15 nm), indicating that such 
magnetism in TI is not from ferromagnetic Eu ions in- 
terdiffused into Sb 2 Te 3 , but from proximity exchange ef¬ 
fect. The free of interdiffusion is also consistent with our 
Transmission Electron Microscopy (TEM) result in Fig. 
§d), where a sharp interface between epitaxially grown 
EuS and Sb 1 . 9 V 0 . 1 Te 3 is developed. 

The magnetization at the interface in proximity 
structures is greatly enhanced when V-dopant is at 
present, from 1.2/is/u.c (Fig. §>. without V-doping) to 
2.7/is/u.c (Fig. i, V-doped), i.e. almost tripled. In 
both cases, the penetration depth of proximity is ~ 1 nm, 
consistent with Bi 2 Se 3 / EuS interface [ 22 ]. Besides, 
the in-plane magnetization of EuS drops dramatically 
near the interface, from ~ 3/ie/u.c without V-dopants 
to ^ 0.5/iB/u.c with V-doping, at fi^H = 5mT. This is 
due purely to magnetic effect instead of interfacial rough¬ 
ness since the ASLD is flat near the TI interface. On the 
contrary, magnetization drop at EuS / AI 2 O 3 interface 
(^ 23 nm) is due to the Stranski-Krastanov growth m, 
leading to a thickness variation and formation of island. 
This is directly confirmed from Z-contrast high-angle an¬ 
gular dark field (HAADF) TEM image (Fig. [5Jd)). At 
higher field fi^H = 0.7 T, an increase of the in-plane EuS 
magnetism is accompanied with a drop of proximity ef¬ 
fect into TI. Since only the perpendicular direction mag¬ 
netism will contribute to the proximity effect pQ, a high 
in-plane guide field tends to align the EuS moment in¬ 
plane and reduce the proximity. 

To understand the origin of the in-plane EuS mag¬ 
netism drop near the interface, we examined the exchange 
bias (EB) of the magnetic hysteresis measurements. Fig. 
[4] plots the results of low-field in-plane hysteresis mea¬ 


surements of a 2 nm EuS / 10 QL Bii.9Vo.lTe 3 hybrid 
heterostructure instead of Sb 2 Te 3 since both belong to 
Bi 2 Se 3 TI family and share very similar crystalline struc¬ 
ture. Despite Sb 2 Te 3 is more suitable for PNR studies 
due to less interstitial V-defects, BUTe 3 is more suitable 
for SQUID due to higher diamagnetic susceptibility. Fig. 
|4ja) shows that the EB can be switched from negative to 
positive by a field of /io H = =t=lT. Fig. [4^b, c) show the 
characteristics of the EB as a function of the resetting 
field, at 5K and 7K, respectively. We adopt the tra¬ 
ditional approach for exchange bias (EB) measurement 
[5211541 at various reset fields, where the EB was initially 
set negative by applying a field of —IT, which was fol¬ 
lowed by a positive resetting field and then measuring 
the low-field hysteresis [35]. This was repeated for re¬ 
setting fields from 0 to + 0.8 T, where the bias is seen 
to shift from = -5 to +60e. Results of the ex¬ 

change biasing strongly suggests the existence of an AF 
structure and a likely magnetic frustration [36] . and is 
quite striking since our system is only composed by two 
strong FMs. The possible configuration derived from EB 
is illustrated in Fig. |4jd), where V-doped TI keeps a per¬ 
pendicular anisotropy, but an interfacial AF structure is 
created to cause the EB. 

To further understand the implication of the results 
in Fig. [4] we develop a phenomenological energy model 
to describe the coupling between the FM and AF. The 
anisotropic energy for bulk EuS can be written as m 

E a n = KiM s t (a\a\^-a\a\ + a\a^)-\-i<i 2 M s ta\a\a\ ( 1 ) 

where oli is the directional cosine along \ th direction, M s 
is the saturation magnetism per area, t is the thickness 
of FM layer and the anisotropic constants k,i = —19.6 Oe 
and k 2 = —4.6 0e at T = 1.3 K [26] . Since our interest 
is in thin film structures with a single symmetry axis ( y- 
axis in Fig. [I] a), eq. 0 could be rewritten as using a 
simplified model for hexagonal and cubic lattice m , 

E an = KiM s tsm 2 (6) + K 2 M s t sin 4 (6) (2) 

where K\ — = —19.6 0e , K 2 = — 2 

1 (}.(} Oe. 6 is the angle between the magnetization and 
the symmetry axis. Since K\ < 0, 6 = n/2 corre¬ 
sponds in the present case for EuS showing easy-plane 
anisotropy within xz- plane. For a thin film, we further 
define Ah,e// = K s /t + K\. We require the surface 
anisotropy constant K s > 0 , since for thinner sample 
Ab, e // will approach zero from the negative side, indi¬ 
cating a rotation of the in-plane easy plane to an out-of¬ 
plane direction, resulting in a magnetic canting which is 
reasonable for compensated thin film interfaces [551 . 

Taking into account the external magnetic field H and 
the FM/AF coupling J, the total energy could be written 
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FIG. 3. (color). PNR fitting profiles of doping-only sample 
Sb 1 . 9 V 0 . 1 Te 3 (a), proximity only sample EuS / Sb 2 Te 3 (b) 
and hybrid heterostructure sample EuS / Sb 1 . 9 V 0 . 1 Te 3 (c). 
The NSLD, MSLD and ASLD denote the nuclear, magnetic 
and absorption scattering length density, respectively. The 
proximity effect appears as finite magnetization signal (blue 
curves) in the region of TI near the interface (~ 15 nm). The 
absorption free feature in this region excludes the possible 
contribution of interdiffused Eu ions. We see clearly that with 
V-doping, the proximity magnetism is enhanced as a bump in 
(c), accompanied with a further suppression of magnetism of 
interfacial EuS (15 — 18nm). (d)HAADF TEM image of the 
EuS / Sb1.9V0.1Te3 hybrid heterostructure. A sharp interface 
between the TI / MI is developed, indicating an epitaxial 
growth of EuS, consistent with (c) for uniformly distributed 
ALSD of EuS. The island-like crystalline facets between EuS 
and AI2O3 cap explains the roughness in (c). 


as 


E = — HM s t sin(0) — J sin(0) 

+ (k! + M s tsm\e) + K 2 M s tsm 4 (6) (3) 

At saturation field configuration 0 = ±7t/2 and consid¬ 
ering the energy extreme, we obtain the bias field and 
coercivity 

/ 2 K 

Hbias = ~Mi ’ H c = 2K 1+ ±K 2 + -^ (4) 

respectively. 

The anisotropic coefficients strongly depend on tem¬ 
perature m In the mean-field approximation, the tem¬ 
perature dependence of anisotropy can be expressed us¬ 
ing the Callen and Callen theory as [38] 

n(n + l) 

K(T) = K( 0)(l-|) 4 (5) 

where n is the order of anisotropy constant, n(K\) = 2 
and n(K 2 ) = 4. Assuming that the Curie temperature of 




+ f IHiS 



FIG. 4. (color). Magnetic measurements of a 2 nm EuS /10 
QL Bi 1 . 9 V 0 .lTe 3 hybrid heterostructure using a Quantum 
Design SQUID magnetometer. (a)In-plane hysteresis showing 
that a negative EB following a set field of - IT, can be switched 
to positive bias by applying a positive resetting field. Inset 
is the out-of-plane magnetic hysteresis of the same sample, 
showing a finite remanent moment, (b, c) EB and coercive 
field as a function of the in-plane resetting field at 5K and 
7K, respectively, (d) Illustration of EB where the moments 
at the interface of EuS are pinned by the exchange-coupled 
moments in the V-doped TI. 


TABLE I. Temperature dependence of the anisotropy con¬ 
stants. Green, Blue, red and black colored values are from 
m , eq. ©> measurements in Fig [4] and eq. @> respectively. 


II 

-4 

1.3 K 

5 K 

7 K 

A'i(Oe) 

-19.6 

-13.10 

-9.96 

K 2 ( Oe) 

+16.6 

7.41 

4.29 

H c ( Oe) 

N/A 

12.7 

4.4 

K s /t( Oe) 

13.8 

9.3 

7.2; 7.1 


EuS is T c = 17 K, we obtain the temperature dependence 
of anisotropy constants as shown in Table [I] One remark¬ 
able feature for this model is that the surface anisotropy 
K s /t calculated from experimental values and eq. © 
coincides with the independent check using eq. ©, giv- 
ing 7.2 Oe vs 7.1 Oe at 7K. Finally, this yields a sur¬ 
face anisotropy K s = 0.0014 erg • cm -2 by assuming a 
2.5 x 10 -5 emu saturation and 5 mm 2 sample area. This 
term is the origin of magnetic canting of interfacial EuS. 

Contrary to the strong T-dependence of anisotropy, the 
bias field i7^ as thus AF/FM coupling constant J has a 
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weak dependence with temperature, indicating an origin 
of FM/AF coupling different from magnetic crystalline 
anisotropym such as the prominent role of Spin-Orbit 
interaction and spin-momentum locking at the surface of 
the TL 

To summarize, we have reported a large enhancement 
of proximity exchange coupling strength in Ml/magnetic- 
TI hybrid heterostructure. This overcomes the major 
disadvantage in MI/TI heterostructures where the prox¬ 
imity coupling strength is considered weak [23]. To our 
knowledge, this is also the first report combining TM 
doping and proximity effect to magnetize TIs. Here, the 
magnetic dopants magnetize the TI surface states, which 
are further coupled to the MI; whereas the magnetic TI 
has a strong perpendicular anisotropy that compensates 
the weakness of MI with in-plane anisotropy. The re¬ 
duction of interfacial magnetism is consistent with the 
exchange bias result, where an AFM-like structure is ar¬ 
tificially created, where a K s = 0.0014 erg • cm -2 surface 
anisotropy is extracted. Despite this value being small 
compared to the stronger examples such as the Au/Co 
interface [39] . this approach provides fruitful insights to 
tailor new magnetic structure at TI/MI interfaces. 
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